Distribution of Singularities in the Cosmic Microwave Background Polarization 
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The polarization of the cosmic microwave background radiation will have a distribution of singu- 
larities and anti-singularities, points where the polarization vanishes for topological reasons. The 
statistics of polarization singularities provides a non-trivial scheme to analyze the polarization maps 
that is distinct from the usual two-point correlation functions. Here we characterize the statistics 
of the singularity distribution in simulated polarization maps, and make predictions that can be 
compared with ongoing and upcoming observations. We use three different characterizations: the 
nearest neighbor distance between singularities, the critical exponent v that captures the scaling 
of total charge q within a closed curve of length L (q oc L"), and the angular two point angular 
correlation functions for singularities of similar and opposite charge. In general, we find that the 
distribution of singularities is random except on scales less than about 10°, where singularities of 
the same charge repel and those of opposite charge attract. These conclusions appear to be ex- 
tremely robust with respect to variations in the underlying cosmological model and the presence of 
non-Gaussianity; the only exception we found are cases where statistical isotropy is grossly violated. 
This suggests that, within the assumption of statistical isotropy, the distribution is a robust feature 
of the last scattering surface and potentially may be used as a tool to discriminate effects that occur 
during photon propagation from the last scattering surface to the present epoch. 



The cosmic microwave background (CMB) 
anisotropics are linearly polarized at the 10% level. 
The polarization, predicted almost four decades ago 
0, has recently been observed @, Q and efforts are 
underway to map it on increasingly smaller scales. The 
polarization is most easily described in terms of so-called 
E and B modes 0, The two-point correlation 

function of E and B maps is typically taken as a statistic 
of choice to represent the polarization properties of 
the map. The two-point function fully describes the 
map if it is Gaussian random and isotropic, otherwise 
higher-order correlation functions are necessary for the 
full description (for a review of CMB polarization, see 
e.g. Ref. 0). 

In this paper we explore another, independent and 
largely unexplored signature of CMB polarization: the 
distribution of singularities 01113. The CMB polariza- 
tion map (denoted P) corresponds to a map from the sky 
- the two dimensional sphere, S 2 - to the space of head- 
less vectors, S 2 /Z 2 , given by the plane and amplitude of 
polarization: 



P : S 2 



S 2 /Z 2 



(1) 



Such maps are known to contain topological features that 
are characterized by points with vanishing polarization, 
known as "singularities" or "defects", and each singu- 
larity carries a topological charge. The total topological 
charge within a closed contour on the sky can be calcu- 
lated as an integral over the contour. This is very sim- 
ilar to Gauss' law that is used to determine the electric 
charge within a closed surface by integrating the elec- 
tric field over the surface. The polarization singularities 
have fundamental charge +1/2. The anti-singularities 
have fundamental charge —1/2. These fundamental sin- 
gularities are shown in Fig. 2] They can be combined 
to form three kinds of double singularities as shown in 



FIG. 1: Fundamental singularities of charge ±1/2. Each 
dash represents the linear polarization of the CMB radiation 
at that point. The x marks the position of the singularity 
where the linear polarization vanishes for topological reasons. 



Fig. [2 "knots" and "foci" which have charge +1, and 
"saddles" with charge —1 7). The total charge in a given 
map is zero for all practical purposes, as discussed in the 
next Section. 

Polarization of the CMB, like the temperature, pro- 
vides an extremely important window to the processes in 
the early universe. In fact, the polarization has some 
advantages over the temperature, and in particular it 
offers a more direct probe of the recombination era at 
large angular scales [lfij . It is therefore well worthwhile 
to consider alternative analyses of the polarization, and 
this motivates us to explore the distribution of singular- 
ities. Our predictions can be tested as soon as polariza- 
tion maps become available. Furthermore, since all of our 
statistics are computed in real (and not Fourier) space, 
sky cuts due to galaxy and other foreground contamina- 
tion will be relatively easy to take into account. 

Polarization singularities have first been described in 
Refs. |7j and |8|. The expected distribution has been de- 
scribed in Ref. [3] based on earlier experience with con- 
densed matter systems. We extend previous work and 
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distance to the nearest neighbor of a singularity. Finally, 
motivated by the analyses of the galaxy distribution on 
the sky, we compute the angular two-point correlation 
function of singularities in the polarization map. This 
work complements the pioneering theoretical discussions 
in Refs. 7, 8, 9] and helps establish the distribution of sin- 
gularities in the CMB polarization as a significant probe 
of the universe. 
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FIG. 2: Singularities of charge ±1, constructed by combining 
the two fundamental singularities shown in Fig. Top row 
shows examples of a "knot" and a "focus", while a "saddle" 
is shown at the bottom, using the nomenclature from Ref. 0. 



In the next section, we describe how we create the 
mock polarization maps. In Sec. IIIII we find the critical 
exponents and the angular correlation function, as well 
as the distribution of distance to the nearest neighbor 
for a vanilla mock map. We investigate the effects of 
non-random phases and non-Gaussianity in Sec. II VI and 
summarize our results in Sec. In the Appendix we 
describe the algorithm for determining the presence of a 
singularity. 



use a highly quantitative analysis of the singularities, de- 
scribe their statistic as a function of pixelization scale 
and the underlying cosmological model, and assess the 
statistical errors for each measurement. 



I. THE DEFINITION OF SINGULARITIES 



II. SCHEME TO PRODUCE MOCK 
POLARIZATION MAPS 



The polarization patterns with fundamental (charge 
±1/2) singularities are shown in Fig.^ The polarization 
map is described by the Stokes parameters Q and U. The 
properties of the polarization under rotations imply that: 

Q = /cos(2a), U = Isin(2a) (2) 

where / is the radiation intensity and a is the polariza- 
tion angle. The ±1/2 singularities are locations around 
which a changes by ±7r. Therefore, given a polarization 
map, we can find the change in a as we go around a 
small closed path and this will tell us whether there is a 
singularity within that closed path. By going around all 
possible paths, we can find all the polarization singulari- 
ties. In the continuum, there would be an infinite number 
of small paths. But in practice, the map is pixelized and 
the change in a is found around circuits defined by neigh- 
boring pixels. The algorithm for finding the singularities 
in a pixelated map is well-known and is described in the 
Appendix. 

In this paper we will be concerned with the distribu- 
tion of polarization singularities. We will characterize the 
distribution in three distinct ways. The first method is 
to find the total charge q within a closed path of length 
L. The mean charge will, of course, be zero because 
one can have positive and negative charges with equal 
probability . Next, we compute the distribution of the 



1 Since the sky is a two sphere, the net charge is +2. This is simply 



To produce mock polarization maps, we proceed as 
follows. We first produce the angular power spectra 
(temperature and polarization) of the CMB using the 
CMBFAST package jn.] ■ Our fiducial model is the stan- 
dard ACDM cosmology with a flat universe with mat- 
ter energy density relative to critical — 0.3, dark 
energy equation of state w = — 1, scalar spectral index 
n s = 1.0, physical matter and baryon energy densities of 
11m/i 2 = 0.127 and fish 2 = 0.021 respectively, and no 
tensor modes (we explore the variations to this model in 
Sec. IIV0 . We can obtain an arbitrary number of statis- 
tically independent maps by generating different sets of 
coefficients af m and af m , consistent with the same under- 
lying power spectra Cf and Cf respectively, and gener- 
ating the polarization map in Healpix |l2j for each set 

separately. Recall that the coefficients a[f n B ^ fully de- 
scribe a given map and, for a Gaussian random map, 
come from a Gaussian distribution of variance Cg E ' B \ 
In some cases we will want to produce polarization maps 
that do not come from Gaussian random ae rn , and in 
those cases we simply input the desired non-Gaussian 

4m B) directly 



due to the topology of the two sphere - it is impossible to comb 
the hair on a two sphere without singularities. This charge is a 
small number compared to the charges that are present due to 
the statistical nature of the polarization map and can be ignored 
for practical purposes. 
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Polarization map for l_max=100, LCDM model Defects in the map 




FIG. 3: Left panel: An example of the mock polarization map. For visual clarity, this map was produced at the resolution 
NSIDE= 256, has power out to ^ max = 50 and contains about 6000 singularities. Right panel: The distribution of singularities 
in the same map. 



III. STATISTICS OF THE SINGULARITY 
DISTRIBUTION 

A. Number of singularities 

We first compute the total number of singularities (or 
charges) N s - ln , positive plus negative, in a given mock 
map. As expected, the number of charges increases with 
the map resolution, just as, for example, the number of 
temperature hot and cold spots increases with resolu- 
tion. N s i n increases roughly as the square of the max- 
imum resolution of the map ^ max so that, for example, 
a map (consistent with our fiducial ACDM model) with 
f max = 100 has about 6000 singularities, while a map 
with £ max = 500 has about 150,000 singularities. Note 
that, even if we are able to measure the polarization with 
infinite resolution, N s i n does not increase indefinitely, but 
levels off when the power in polarization becomes negli- 
gible. For a standard ACDM model, the EE power spec- 
trum has power all the way to I ~ 2000. However, this 
implies that the total number of singularities is likely to 
be larger than a million, making the analysis of such a 
map challenging. Fortunately, we do not need to worry 
about this issue, or wait for polarization experiments that 
will reach scales this small, such as PLANCK, in order to 
explore the distribution of singularities: the tests we pro- 
pose can be performed for polarization maps covering any 
range of angular scales, and statistics from the measured 
map of any given resolution can be compared to Monte 
Carlo tests with mock maps of the same resolution. 



B. Scaling of RMS charge 

We would now like to quantify how the variance in the 
number of singularities increases with the area covered on 



any given map. The root-mean-square (RMS) fluctuation 
of the charge within a closed path of length L is expected 
to be: 

<j(q(L)) ^ ((q(L) - q(L)f) 1/2 = aL» (3) 

where q(L) is the total (positive plus negative) charge 
within L, q(L) is its mean among the different paths of 
the same length, a is a system-dependent constant and 
the critical exponent v is expected to be 0.5 J9|- Using 
our numerical analysis of mock data we will be able to 
predict both a and v together with error bars. 

To compute the RMS of charge per ring, we create 
a polarization map, and pick a point on it that we call 
the North Pole. In the Healpix representation of the 
map, each ring of angle 9 from the North Pole has length 
L = 2nsmQ, and we find q(L) for all rings using the 
procedure described in the Appendix. We then rotate the 
north pole of the map in a random direction (i.e. assign 
it to a new, randomly chosen point on the sphere) and 
repeat the computation of q(L) for each ring. We repeat 
this procedure hundred times in order to obtain sufficient 
statistics and compute a(q(L)). The mean charge per 
ring, averaged over all rotations, is nearly zero, while the 
fluctuations around the mean is what we are interested 
in. 

Figure0](top panel) shows the scaling of a(q(L)) with 
L for maps of two resolutions £ max = 100 and 500, cor- 
responding to polarization having power down to scales 
of ~ 1° and ~ 0.4° respectively. In both cases we have 
fixed the pixelization of the map to the Healpix parame- 
ter NSIDE= 256, corresponding to pixels of about 0.25° 
on a side. In these and many other cases we have ex- 
plored, the RMS of charge scales as L v where the criti- 
cal exponent v is nearly 0.5, in agreement with expecta- 
tions. Furthermore, the exponent is independent of the 
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FIG. 4: Top panel: Root-mean-square of the total charge per 
ring as a function of the ring's length L. The two datasets 
correspond to polarization maps with resolution ^ max = 100 
(bottom) and ^ max = 500 (top), while the dashed curves de- 
note the linear fit in the log-log coordinates. The distributions 
are consistent with a(q(L)) oc V with 0.45 < v < 0.55. Bot- 
tom panel: the dependence of v on the maximum resolution 
of the map ^ ma x (^max is the maximum multipole where the 
map has power). 
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FIG. 5: Histogram of the distance from any given singularity 
to the nearest neighbor of the same charge (black-solid) and 
opposite charge (red-dashed). Note the effects of repulsion 
between charges of the same sign, and attraction of charges 
of the opposite sign. 



C. Distance to the nearest neighbor 

Another statistic that we explore is the distance to the 
nearest neighbor from any given singularity. The his- 
togram of the distances to the nearest neighbor is shown 
in Fig. [SJ where for computational convenience we have 
assumed a map with £ max — 100 which has a total of 
about 6000 singularities. The average distance to the 
nearest neighbor is slightly above 2° and can be roughly 
predicted from the total number of singularities (~ 6000 
singularities in ~ 40000 degrees on the sky). However, 
the histograms of the nearest neighbor of the same charge 
and that of the opposite charge are different, and show 
that opposite charges attract and similar charges repel. 
For example, the mean distance to the nearest neighbor 
of the same charge is (3.24 ±0.01)°, while distance to the 
neighbor of the opposite charge is (2.69 ± 0.01)°. 



map pixelization NSIDE 2 . However, we find that the ex- 
ponent slightly decreases with map resolution, moving 
from » 0.55 for € max = 100 to » 0.45 for £ max = 500; see 
the bottom panel of Fig. 0] In all cases the error in the 
exponent is about 0.01. 



2 It is important to use NSIDE large enough to capture the res- 
olution of the map and avoid pixel effects. This corresponds to 



NSIDE> . 



,/4. 



D. Angular power spectrum of singularities 

The third characterization of the singularity distribu- 
tion we suggest and explore is the two point angular cor- 
relation function w(9). The angular two point function is 
simply the excess probability, on top of expectation due 
to random distribution, of finding one singularity at an 
angular location 0\ and another one at location 0i 



dP(0 1 ,0 2 ) = dP(e 1 )dP(d 2 )(i + w(e 1 -e 2 )) (4) 
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dP(6 1 )dP(e 2 )(l + w(6)) 



(5) 



where 8 = \ 6\ — 62 | and the second line assumes sta- 
tistical isotropy. In addition to the angular correlation 
function for all singularities, we can find individual w{9) 
for singularities of positive (or negative) charge in a sim- 
ilar way. 

The angular correlation function is one of the standard 
tools to describe the angular clustering of galaxies, and 
has been thoroughly explored and used during the past 
three decades. Computing w(9) for the galaxy distribu- 
tion, however, typically involves various practical prob- 
lems, most important of which is the "selection function" 
of the survey, having to do with magnitude cuts and im- 
perfect coverage across the field of view. The application 
we are considering here is vastly simpler, since the singu- 
larities are discrete, well-defined and easily computable 
features, and unlike the galaxies they are located at the 
same radial distance. Furthermore, we are simulating 
the polarization pattern on full skies and do not need to 
worry about edge-effects due to incomplete sky coverage. 
Therefore, a simple estimator for w(9) will suffice. We 
adopt the Peebles-Hauser estimator 



w(6) 



AUap \RR(6) 



(6) 



where DD(9) is the number of singularity pairs separated 
by a distance larger than 9 — d9 but smaller than 9 + d9, 
RR(9) is the number of pairs from a random (Poisson- 
distributed) map being in the same distance interval, and 
-'Vmap and iV ran d are the total numbers of singularities in 
the two maps. In other words, the angular correlation 
function measures the excess clustering over that pre- 
dicted by random distribution on any given scale. Here 
we use d9 = 5° , which is larger than our pixelization scale 
and thus avoids any edge-effects due to pixelization. Wc 
have tried several other values of d9 and obtained consis- 
tent results. 

Figure shows the angular two-point correlation func- 
tion for the map with £ max — 100. In each case we com- 
pute w(9) for 10 statistically independent maps, and plot 
the mean and standard deviation of measurements at all 
angular scales. Note that this procedure correctly ac- 
counts for the increase in the error bars at large angular 
separations due to sample (or cosmic) variance. 

We find that w(9) is remarkably consistent with zero 
at all angular scales greater than about 10°! This implies 
that the singularities are Poisson distributed over most 
angular scales. Also note that the singularity-singularity 
w(9) is negative at scales < 10°, while the singularity- 
antisingularity w(9) is positive over the same range. This 
illustrates the fact we discussed in Sec. lIII CI that charges 
of the same sign repel and those of opposite sign at- 
tract. Finally, w(9) for all singularities (irrespective of 
their charge) shows an overall decrease at 9 < 10°. We 
have repeated this analysis by varying the maximum res- 
olution of the map £ maK and found consistent results: the 
behavior of w{9) is qualitatively similar, while the error 
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FIG. 6: Angular correlation function of the distribution of 
singularities, showing the excess probability of clustering over 
that predicted by the random distribution. The overall dis- 
tribution of singularities is consistent with random, except 
at small scales (< 10°) where singularities of opposite charge 
attract while those of equal charge repel. 



bars decrease with increasing £ n 
the number of singularities. 
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IV. EXPLORING DEPARTURES FROM THE 
STANDARD COSMOLOGICAL MODEL 

It is important to determine how the singularity distri- 
bution depends on the physical input such as the primor- 
dial fluctuation spectrum, the cosmological parameters, 
or Gaussianity of cosmological seed perturbations. In 
particular, previous work 0, @ emphasized that the dis- 
tribution and type of singularities may be a promising 
way of probing the Gaussianity of initial conditions. 

With this in mind we have created mock maps of the 
CMB polarization using several different cosmological 
models and characterized the singularity distribution in 
each. In particular, we have tried several extreme possi- 
bilities, some of which are already ruled out by cosmo- 
logical observations: 

• Maps based on primordial power spectrum with 
very significant tensor modes (the ratio of tensor 
to scalar perturbations at the CMB temperature 
quadrupole of T/S = 10). 

• Maps with a strongly tilted primordial power spec- 
trum with either less or more power on small scales: 
we alternatively assumed a scalar spectral index 
n s = or n s = 2. 
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• Maps that are strongly non-Gaussian: we assumed 
the real and imaginary parts of the coefficients 

a \m B ^ t° be sampled from an exponential distribu- 
tion while keeping zero mean and variance equal to 
/2. In other words, we have adopted a highly 
skewed distribution of the coefficients a^ B ^ . 

Remarkably we find that all statistics we considered 
are largely unaffected. The critical exponent is un- 
changed within the errors, and so is its dependence on 
f mffl . The two point correlation function w(9) is also 
unchanged, being consistent with zero except on scales 
less than 10°. The total number of singularities does 
change for the alternative cosmological models, which is 
to be expected since the total power in the map is af- 
fected each time. However, statistics of the distribution 
of singularities are unchanged. 

In particular, we find it very interesting that the re- 
sults are insensitive to the model for non-Gaussianity we 
assumed. We have explored this further for several other 
classes of variations to the standard Gaussian/isotropic 
assumption, and found deviations from the results de- 
scribed in Sec. IIIII only in cases where the statistics of the 
map were modified enough that the statistical isotropy 
was grossly violated. For example, maps where the only 
nonzero a^ m ' coefficients were those with m = showed 
deviation from results in Sec. IIIII however, inspection of 
polarization maps in these cases indicate that such mod- 
ifications lead to huge violations of isotropy in the map 
that are easily detectable with almost any reasonable sta- 
tistical test. Conversely, more subtle modifications of 
the power spectra (for example, setting the quadrupole 
of temperature and polarization to zero) produced the 
same results as our fiducial model. Therefore the charac- 
teristics of the singularity distribution seem to be robust 
features that are insensitive to the cosmological inputs at 
the last scattering surface. 

The robustness of the singularity distribution holds ad- 
vantages as well as disadvantages. One disadvantage is 
that by examining the distribution of the singularities in 
the actual data, it is unlikely that we will uncover some- 
thing about the physics of the primordial fluctuations 
(unless statistical isotropy is violated). The advantage 
is that since the distribution is robust, any distortions 
in it must come during propagation of the photons from 
the last scattering surface to us. Hence the distribution 
can be used as a probe of cosmology at redshifts smaller 
than 1000. For example, weak gravitational lensing could 
distort the Poissonian nature of the angular correlation 
functions; however, this effect will operate only at small 
angular scales, as lensing distortions are typically a few 
arcminutes. 



V. CONCLUSIONS 

We have explored the statistics of the distribution of 
singularities in the CMB polarization maps. The exis- 



tence of singularities are a generic prediction for a head- 
less vector field on a sphere. Their distribution, however, 
has not been explored in the past except for some generic 
scaling arguments. Here we have provided a quantitative 
analysis of the distribution of singularities with positive 
and negative charge, and argued that the singularities 
provide an additional probe of the conditions at last scat- 
tering, largely independent of the usual two-point corre- 
lation functions of temperature. 

We found that the singularities are distributed ran- 
domly at scales > 10°. The angular two-point correlation 
function of singularities vanishes at those scales, while 
the total charge within a closed path of length L scales 
as L v with v ~ 0.5. On scales < 10°, however, charges 
of the same sign repel while those of opposite charge at- 
tract. The attraction and repulsion are manifested both 
in the angular two-point correlation function and in the 
distribution of the nearest-neighbor distance. 

Perhaps surprisingly, we found that the aforemen- 
tioned results are very robust with respect to variations 
in the underlying cosmological model assumed to create 
the mock maps. Changes in the tensor to scalar ratio, 
scalar spectral index, and non-Gaussian distributions of 
the a; m 's all leave the singularity distribution unchanged. 
Only violations of statistical isotropy affected the distri- 
bution of singularities. It is still possible that there are 
some other forms of non-Gaussianity that can affect the 
singularity distribution that we have not explored. While 
this is impossible to rule out, it seems more likely to us 
that the distribution of singularities at the last scatter- 
ing surface is described precisely by the Poissonian dis- 
tribution and other characteristics we have found. This 
implies that any observed deviations of the actual map 
from these distributions will have to be due to line of 
sight effects. Most importantly, gravitational lensing of 
the CMB by the large-scale structure can cause changes 
in the singularity distribution. However, the lensing op- 
erates mostly on small scales, with typical deflections of a 
few arcminutes and coherence of < 10°. Therefore, lens- 
ing may affect the statistics of the singularity distribution 
only at small scales. We hope to explore this signature 
in future work. 
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FIG. 7: To calculate the topological winding along a contour 
r - here shown as a hexagon - we follow the phase a, always 
choosing the shortest path on the a circle, ai denotes the 
value of the phase at point i on F. The net change in a/2ir 
as r is traversed is the winding. In the figure, the change in 
a is — 7r since we go around the a circle once in the counter- 
clockwise direction. Note that the full circle in a corresponds 
to an angle of only 7T, not 2n. Hence the topological charge 
within F in the drawn example is —1/2. 



APPENDIX A: CALCULATION OF THE 
WINDING NUMBER 

Here we will describe how the winding of polarization 
around a contour T, i.e. net topological charge within T, 
is calculated. 

The polarization map specifies Q{6,(f>) and U{8,4>), 
where and 4> are galactic latitude and longitude. From 
Q and U we can determine the polarization angle a: 

a(e,<P) = ^- 1 (j^j (Al) 

Now we want to know the change in a as we go around 
the closed contour T. 



Any map of the CMB polarization will be pixelized. 
Hence only an average value of a within each pixel will 
be available to us. Then the change in a in going from 
pixel i to a neighboring pixel i + 1 is given by: 

Sat = —on + p (A2) 

where (3 is defined as follows: 

(3 = 0, if \a i+1 -c^l <tt/2 (A3) 

f3 = +7r , if a l+ i - a. L < -tt/2 (A4) 

(3 = — 7T , if a l+1 - on > +7r/2 (A5) 

In other words 5a is the shortest path from ai to a^+i 
around the circle defined by a 6 [— 7r/2, +tt/2] (see 
Fig. EJ). Note that a = —tt/2 yields the same Q and U 
'as a — +tt/2 and hence these two points are identified. 
The winding, Aa, around the contour T is now simply: 

Aa = 5a i ( A6 ) 

i 

where the sum is over all the discretized steps that define 
r. The resulting topological charge is then Aa/(27r). 

This scheme to find the windings has one small mod- 
ification since the sky is S 2 . In this case, since the po- 
larization angle a is defined with respect to the lines of 
latitude and this definition breaks down at the North 
and South poles, we need to modify the scheme for any 
contour that contains the North (or the South) pole. In 
that case, 1 must be added (per enclosed pole) to the net 
topological charge within the contour. This also ensures 
that the total topological charge on the sky is +2 9]. 
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